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ABSTRACT 

We report on 10 yr of monitoring of the 8.7-s Anomalous X-ray Pulsar 4U 0142+61 using 
the Rossi X-Ray Timing Explorer (RXTE). This pulsar exhibited stable rotation from 2000 until 
February 2006: the RMS phase residual for a spin-down model which includes i/, 0, and i> is 2.3%. 
We report a possible phase-coherent timing solution valid over a 10-yr span extending back to 
March 1996. A glitch may have occured between 1998 and 2000, but it is not required by the 
existing data. We also report that the source's pulse profile has been evolving in the past 6 years, 
such that the dip of emission between its two peaks has been getting shallower since 2000, almost 
as if the profile is recovering to its pre-2000 morphology, in which there was no clear distinction 
between the peaks. These profile variations are seen in the 2—4 kcV band but not in 6—8 keV. 
Finally, we present the pulsed flux time series of the source in 2—10 keV. There is evidence of a 
slow but steady increase in the source's pulsed flux since 2000. The pulsed flux variability and 
the narrow-band pulse profile changes present interesting challenges to aspects of the magnetar 
model. 

Subject headings: pulsars: individual(4U 0142+61) — stars: neutron — X-rays: stars — 



1. Introduction 

1.1. Anomalous X-ray Pulsars (AXPs) 

The existence of magnetars - young, isolated 
neutron stars powered by the decay of an ultra- 
high magnetic field - is now well supported by sev- 
eral lines of evidence (Woods & Thompson 2006). 
There are at least two fiavors of magnetars: SGRs 
and AXPs. They both exhibit: X-ray pulsations 
with a luminosity of iQ-^^'^^e ^^.g g-i^ periods rang- 
ing from 5-12 s, period derivatives of 10""'^^ — 10~^^, 
and surface dipolar magnetic fields of 0.6—7 x 
10^* G. In the magnetar model, the pulsed X- 
rays are the result of a combination of surface 
thermal emission and resonant scattering in the 
magnetosphere (Thompson, Lyutikov & Kulkarni 
2002). For more on AXPs and SGRs, see reviews 
by V. Kaspi and S. Mereghetti (this volume). 
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1.2. AXP 4U 0142+61 

4U 0142+61 is an 8.7-s AXP. It has P = 0.2 x 10"" , 
implying a surface dipole magnetic field of 1.3 x 10^^ G^. 
It is known to pulsate in the optical band (Kern & 
Martin 2002; Dhillon et al. 2005), and has been de- 
tected in the near-IR (HuUeman, van Kerkwijk & 
Kulkarni 2004), in the far-IR (Wang, Chakrabarty 
& Kaplan 2006), and in hard X-rays (Den Hartog 
et al. 2006; Kuiper et al. 2006). It has a soft X-ray 
spectrum well fitted by a combination of a black- 
body and a power law (see, for example. White 
et al. 1996). AXP 4U 0142+61 rotates with high 
stability (Gavriil & Kaspi 2002). However, Morii, 
Kawai & Shibazaki (2005) reported a timing glitch 
in 1999 on the basis of an ASCA observation in 
which the value of the frequency is marginally 
discrepant with that predicted by the ephemeris 
reported by Gavriil & Kaspi (2002). 

Here we report on continued RXTE monitoring 



^Magnetic fields here are calculated via B = 3.2 X 
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observations of this source. Our observations are 
described in Section 2. Our timing, pulsed mor- 
phology, and pulsed flux analysis are presented, 
respectively, in Sections 3, 4, and 5. 

2. Observations 

We used 136 RXTE/PCA observations of var- 
ious lengths in our analysis: a) 4 very closely 
spaced RXTE Cycle 1 observations, h) lA short 
Cycle 2 observations spanning a period of a year, 
c) 1 Cycle 3 observation (followed by a 2-yr gap 
with no observations), d) 118 observations taken 
regularly from 2000 to 2006 as part of a long-term 
monitoring program spanning RXTE Cycles 5 to 
10. 

For each observation, photon arrival times were 
barycentered and binned with 31.25-ms time res- 
olution. 

3. Phase-coherent Timing 

Each binned time series was folded at the 
pulse period. Resulting pulse profiles were cross- 
correlated with a high S/N template. This re- 
turned an average pulse time of arrival (TOA) for 
each observation corresponding to a fixed pulse 
phase. The pulse phase at any time can be ex- 
pressed as a Taylor expansion polynomial. The 
TOAs were fitted to the polynomial using the pul- 
sar timing software package TEMPO^. 

We report a phase-coherent timing solution 
that spans the post-gap {i.e. after 2000) 6-yr pe- 
riod up until February 2006 (MJD 53787) includ- 
ing all data in RXTE Cycles 5—10. The param- 
eters of our best-fit spin-down model are shown 
in Table 1. The phase residuals are shown in 
Figure 1. The best-fit post-gap ephemeris does 
not, however, fit the pre-gap TOAs well. Figure 2 
shows a clear systematic deviation in the pre-gap 
residuals obtained after subtracting the post-gap 
ephemeris. This could indicate that a glitch oc- 
cured at some time during the gap. However, by 
using 6 frequency derivatives, we found a possible 
ephemeris that fits the entire Cycle 1 to Cycle 10 
range (see Table 1, Figure 3). 

The existence of our overall ephemeris cannot 
rule out the possibility of the glitch having oc- 
cured in 1999 (Morii, Kawai & Shibazaki 2005): 



Year 

2000 2001 2002 2003 2004 2005 2006 




52000 52500 53000 53500 
Modified Julian Days 



Fig. 1. — Arrival time residuals for 4U 0142-1-61 for 

the post-gap period, using the post-gap ephemeris 
given in Table 1. The residuals have RMS 2.3% of 
the pulse period. 



Year 

2000 2002 2004 2006 



51000 52000 53000 

Modified Julian Days 



Fig. 2. — Arrival time residuals for 4U 0142-1-61 
for all Cycles using the post-gap ephemeris. 



^http://pulsar. princeton.edu/tempo. 
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Table 1 

Spin Parameters for 4U 0142+61 





Pre-Gap Ephemeris'' 


Post- Gap Ephemeris 


Overall Ephemeris 


Parameter 


Spanning 


Spanning 


Spanning 




Cycles 1 to 3 


Cycles 5 to 10 


All Cycles 



MJD range 






50170.693-50893.288 


51610.636-53787.372 


50170.693-53787.3 


TOAs 






19 


118 


137 


V (Hz) 






0.115096869(3) 


0.1150969337(3) 


0.1150969304(2) 


V (10"" Hz s" 






-2.659(3) 


-2.6935(9) 


-2.6514(7) 


V (10-2^ Hz s" 








0.417(10) 


-1.7(2) 


d^vjiM:^ (10"^^ 
d^vjdt^ (10-'''' 
d^vjdt'' (10-^« 


Hz s" 








3.62(12) 


Hz 








8.7(3) 


Hz s" 








-5.01(13) 


d^vjdf (IQ-^* 


Hz s" 


-6) 






6.6(4) 


Epoch (MJD) 






51704.000025 


51704.000025 


51704.000000 


RMS residual 






0.019 


0.023 


0.019 



Year 

1996 1996 2000 2002 2004 2006 




51000 52000 53000 

Modified Julian Days 



Fig. 3. — Arrival time residuals for all Cycles using 
the overall ephemeris (see Table 1). The residuals 
have RMS 1.9% of the pulse period. 



if a fully recovered glitch with a short relaxation 
time occured in the RXTE/PCA observing gap 
between Cycles 3 and 5, only a random phase 
jump would be observed. To investigate this, we 
added an arbitrary but constant time jump to all 
the post-gap TOAs. We were still able to find a 
new ciphemeris that connected the TOAs through 
the two-year gap. This indicates that our overall 
ephemeris is not unique. Hence, we cannot rule 
out the possibility of a phase jump between Cy- 
cles 3 and 5, and therefore a glitch in 1999 cannot 
be ruled out. However, it is also not required by 
the data. Indeed, the frequency of the discrepant 
ASCA observation reported by Morii, Kawai & 
Shibazaki (2005) is consistent with the frequency 
predicted by the overall ephemeris shown in Ta- 
ble 1 to within 2a. 

4. Pulse Profile Changes 

To search for pulse profile changes, the phase- 
aligned profiles were averaged for each RXTE Cy- 
cle. This was done in three energy bands. The 
average profiles in all bands are presented in Fig- 
ure 4. The ratios of the Fourier amplitudes of the 
pulse profiles in all bands as a function of RXTE 
Cycle are presented in Figure 5. 

In Figure 4, the pulse profile changes are clear: 
for 2-10 and 2-4 keV, in Cycles 1 and 2, the 
smaller peak is not well defined. After the two- 
year gap, in Cycle 5, the dip between the peaks is 
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Pulse Phase (cycles) Pulse Phase (cycles) Pulse Phase (cycles) 
2-10 keV 2-4 keV 6-8 keV 



Fig. 4. — Average pulse profiles in all RXTE Cy- 
cles in 2—10, 2—4, and 6—8 keV. In a given band, 
the different profile qualities arc due to different 
net exposure times. The Cycle number is shown 
in the top right corner of each plot. 
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Fig. 5. Ratios of the Fourier amplitudes of the 
pulse profiles in all three energy bands as a func- 
tion of RXTE Cycle. Top: ratio of the 2nd ampli- 
tude (A2) to that of the fundamental (Al). Mid- 
dle:(A3/Al). Bottom: (A4/A1). 
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much more pronounced. The peaks start to merge 
back in subsequent Cycles almost as if the profile is 
recovering to its original morphology. In 6—8 keV, 
the smaller peak is not as obvious, indicating that 
it has a softer spectrum relative to the larger peak. 

Given the spectrum of the source (see White 
et al. 1996), which is fitted to a two-component 
model consisting of power law and thermal emis- 
sion, the 2—4 keV band includes both thermal and 
power-law photons while the 6—8 kcV band con- 
tains negligibly few thermal photons. In Figure 5, 
the fact that the ratio of the harmonics is dropping 
only in 2—4 keV suggests that only the thermal 
component of the spectrum is evolving. 

5. Pulsed Flux 

For each observation, the pulsed flux, Frms 
(in counts/s/PCU), was calculated by taking the 
square root of the average of the squares of the 
deviations from the mean number of counts in the 
pulse profile. We omitted PCU from this analy- 
sis because of the uncertainties in its response due 
to the loss of the propane layer. 

The pulsed flux scries for 4U 0142+61 shows a 
slow but steady increase since 2000 in 2—10 keV 
(see Fig. 6). There are hints that the change is also 
present in 2— 4 kcV and not in 6— 8 kcV but our 
statistics do not let us confirm this. We verified 
there are no comparable trends in the long-term 
light curves of the other AXPs observed as part of 
this monitoring program. 

From Thompson, Lyutikov & Kulkarni (2002), 
increases in the twist angle of the field lines in 
the magnetosphere can cause luminosity increases 
as well as pulse profile changes. It is tempting 
to also attribute the increase in the pulsed flux 
of this source to an increase in the twist angle. 
However, Thompson, Lyutikov & Kulkarni (2002) 
also predict that an increase in the twist angle 
is accompanied by an increase in the spin down 
rate, which is not what we are currently observing: 
in the post-gap cphcmeris presented in table 1, i> 
is positive, i.e. the magnitude of the spin down 
rate is decreasing. Thus, another explanation is 
needed. 

Finally, it is interesting to note that in April 2006, 
after the last RXTE Cycle included in this analy- 
sis, the pulsar appears to have entered an extended 
active phase: a single burst accompanied by a 



pulse proflle change was detected from the pulsar 
on Aprfl 06 (Kaspi, Dib & Gavriil 2006). A series 
of four bursts was later detected on June 25 (Dib 
et al. 2006). This interesting turn of events is 
presently under careful study. 
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Fig. 6. — Pulsed flux time series for 4U 0142+61 
in the 2-10 keV band. 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.2. 
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